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ABSTRACT 


A  control  system  was  designed  to  attenuate  vertical 
accelerations  for  the  XR-3  captured  air  bubble  type  surface 
effect  ship  using  linear  regulator  techniques  applied  to  the 
simplified  nonlinear  equations  of  motion,  A  pressure  lift- 
only  model  was  used  to  represent  the  craft  vertical  heave 
motion  and  was  linearized  around  the  steady  state  operating 
point.  Model  validation  was  obtained  through  analysis  of 
the  frequency  spectrum.  State  variable  feedback  was  used  to 
determine  a  set  of  optimal  control  gains  that  would  reduce 
the  magnitude  of  the  heave  acceleration  during  operation 
under  simulated  sea  input  conditions. 
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I .  INTRODUCTION 


An  important  design  consideration  for  the  high  speed 
captured  air  bubble  type  surface  effect  ship  is  the  vertical 
acceleration  associated  with  operation  in  a  random  ocean 
environment  [Ref,  1],  The  vertical  motion,  heave,  is  caused 
by  waves  changing  the  volume  of  the  plenum  producing  pres¬ 
sure  variations  that  result  in  vertical  craft  accelerations. 

Studies  of  conventional  hull  ships  show  that  the  strongest 
vertical  plane  accelerations  are  experienced  at  encounter 
frequencies  close  to  the  ship's  natural  frequency.  Reference 
2  states  that  "The  most  severe  and  objectionable  motions  (wet 
decks,  slamming,  high  accelerations)  result  from  the  responses 
to  the  components  that  are  near  synchronism  with  the  ship's 
natural  pitching  and  heaving  periods".  Encounter  frequency 
is  a  function  of  wave  frequency,  craft  speed  and  direction. 
Since  most  of  the  weight  of  the  surface  effect  ship  is  sup¬ 
ported  by  a  bubble  of  air,  large  accelerations  can  occur  that 
affect  habitability  and  sea  keeping  ability  while,  operating 
in  high  seas . 

Methods  of  attenuating  vertical  motion  have  been  the 
subject  of  past  studies.  A  conventional  design  by  Rohr  Marine 
Inc.  [Ref,  3]  utilized  accelerometer  output  to  monitor  verti¬ 
cal  motion  variations.  These  signals  are  passed  through 
compensation  circuits,  filtered  and  applied  to  actuators 


which  control  air  flow  rates  into  the  plenum.  State  variable 
techniques  are  applied  to  describe  the  system  and  control 
functions.  The  physical  constraints  of  the  actuators  are 
considered  with  regard  to  their  inherent  nonlinearities  and 
limitations,  The  sea  input  is  modeled  using  a  wave  height 
power  spectral  density  procedure. 

Through  the  use  of  a  towing  tank  model,  and  data  gathered 
from  the  100B  test  craft,  an  investigation  was  conducted  to 
measure  the  characteristic  response  of  the  SES  while  operating 
in  an  actual  sea  state  environment  [Ref.  4],  The  model  was 
subjected  to  simulated  sea  state  inputs  of  various  severity  in 
an  effort  to  establish  the  linear  characteristics  of  the 
vertical  motion  and  determine  the  damped  natural  frequency  of 
response.  A  series  of  curves  showing  response  amplitude  for 
both  C.G.  accelerations  and  pitch  angle  versus  encounter  fre¬ 
quency  were  developed  for  a  range  of  speed  and  sea  state 
operation . 

In  Ref.  1,  the  frequency  response  approach  was  used  in  an 
analysis  which  examined  the  effects  of  changing  plenum  air 
flow  rates  on  center  of  gravity  accelerations  using  the  six 
degree  of  freedom  simulation  program  of  the  XR-3  [Ref.  5]. 

This  study  developed  a  linear  simplified  model  in  close  agree¬ 
ment  with  the  6  D.O.F.  model  and  in  addition  used  signifi¬ 
cantly  less  computer  time, 

In  Ref.  6,  Boggio  conducted  a  study  on  the  effect  of  the 
use  of  a  flexible  membrane  in  the  plenum  chamber  itself. 

Using  the  loads  and  motions  program  for  the  XR-3,  the  study 
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reduced  computation  time  through  the  inclusion  of  pressure 
rate  and  volume  rate  equations.  Membrane  insertion  served  to 
effectively  reduce  the  magnitude  of  pressure  variations  in 
the  plenum  thereby  attenuating  vertical  motions. 

Grant  [Ref,.  7]  used  simple  models  for  the  control  devices 
applied  to  the  loads  and  motions  program  of  the  100B  SES 
craft.  From  monitored  plenum  pressure  variations  and  heave 
accelerations,  control  signals  were  developed  by  filter  de¬ 
sign  methods ,  to  vary  the  fan  speed  and  vent  louver  openings 
in  order  to  reduce  the  vertical  accelerations  and  pitch  motion. 

This  thesis  applies  the  linear  regulator  design  technique 
to  determine  a  set  of  optimal  state  feedback  gains  to  control 
the  air  flow  rate.  The  equations  of  motion  used  in  the  design 
were  derived  from  a  simplified  pressure  lift  only  model  lin¬ 
earized  around  the  steady  state  operating  point.  Design  re¬ 
sults  were  evaluated  using  both  the  linear  and  nonlinear 
models  . 
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II.  PRESENTATION  OF  THE  SYSTEM 


Two  major  forces  which  support  the  weight  of  the  surface 
effect  ship  are  the  buoyancy  and  plenum  pressure  forces. 

The  buoyancy  forces  (F^)  are  exerted  on  the  craft  by  dis¬ 
placement  of  water  by  the  underwater  volume  of  the  rigid 
sidewalls.  The  plenum  pressure  lift  force  (Fp)  is  exerted 
on  the  craft  by  the  captured  cushion  of  air.  Flexible  bow 
and  stern  seals  along  with  the  rigid  sidewalls  and  lift  fans 
form  the  captured  air  bubble.  Relative  contributions  to 
craft  vertical  plane  motion  by  Fp  and  Fp  depend  upon  the 
encounter  frequency  (u  ).  At  low  frequencies,  F.  responds 
to  the  rise  and  fall  of  the  sea.  At  high  we ,  Fp  attenuates 
because  incremental  volume  changes  of  the  immersed  hull  are 
reduced.  The  Fp  force  responds  inversly  to  changes  in  plenum 
volume  (Vg)  according  to  the  adiabatic  law  of  pressure 
variations  within  the  plenum  chamber  [Ref.  1],  Volumetric 
changes  in  the  plenum  depend  upon  the  sea  state  and  encounter 
frequency.  As  is  shown  in  Ref.  1,  the  effect  of  the  volume 
changes  produces  an  increase  in  the  Fp  force  as  the  encounter 
frequency  is  increased. 

In  the  system  developed  in  this  study,  the  buoyancy  terms 
have  been  omitted  in  an  effort  to  simplify  the  system  further. 
This  simplified  model  is  referred  to  as  the  "pressure  lift 
only"  system.  Transfer  functions  of  acceleration,  volumetric 


plenum  air  mass  to  input  waveheight  and  system  state  equa¬ 
tions  were  developed  from  the  linearized  form  of  the  pres¬ 
sure  lift  only  system. 


The  basic  assumptions  which  form  the  basis  for  the  pres¬ 
sure  lift  only  model  are  presented  below: 

1.  The  center  of  pressure  (C^)  is  directly  under  the 
center  of  gravity  (C.G,). 

2.  The  sidewalls  are  of  uniform  cross  section  and  are 
symmetrical  about  the  C.G. 

3.  The  effect  of  pitch  moments  on  the  vertical  plane 
forces  are  neglected. 

4.  The  only  variable  lift  forces  are  those  forces 
caused  by  changes  in  plenum  pressure  (P^) . 

5.  The  craft  is  at  a  constant  speed. 

6.  The  air  leakage  rate  from  the  rear  seal  is  constant. 

7.  The  craft  weight  is  supported  by  two  component 
forces  in  the  following  approximate  proportions; 

Sidewalls:  20%  (fixed  force) 

Captured  Air  Bubble:  80%  (variable  force) 

8.  Encounter  frequencies  of  interest  lie  within  the 
range  of  17<oae<100  radians/sec. 


A.  EQUATIONS  OF  MOTION 

The  development  of  the  nonlinear  simplified  model  for 
vertical  motion  is  presented  in  Ref.  1  and  the  block  diagram 
appears  in  Appendix  A,  The  development  of  the  pressure  lift 
only  nonlinear  model  follows  in  a  similar  manner.  The  equa¬ 
tions  of  motion  differ  in  that  the  F^  force  is  considered 


constant . 


The  application  of  the  given  assumptions  to  the  Newtonian 
Laws  for  vertical  motion  of  a  rigid  body  system  leads  to  the 
following  equations; 

1)  Orifice  Leakage  Rate 


2)  Fan  Map  Input  Rate 

q.  =  n[Q .  (o)  -  K  P,  ] 
^in  q  b 


sec 


sec 


(2.A.1) 


(2. A. 2) 


3)  Absolute  Plenum  Pressure 


P  _  P  (  Mb  /  lbs 

pb  ~  a  (  VTF~  }  TTT 

bMa  ft 


4)  Plenum  Air  Flow  Rate 


x  lbs 

^b  ~  pa  ^in  "  ^out  sec 


5)  Heave  Acceleration 


..  T7  F_  Ft  ' 

7  -  ,W  P  _1}  ft 

L  "  VM  M  “  M  sec 


6)  Heave  Velocity 


Z  = 


Sl 


dt 


ft 

sec 


(2. A. 3) 


(2. A. 4) 


(2. A. 5) 


(2, A. 6) 
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7)  Heave  Displacement 


Z  =  J  z  dt 


ft 


8)  Plenum  Pressure  Lift  Force 


F  =  A,  P,  lbs 

p  b  b 


9)  Buoyancy  Lift  Force 


F,  =  C2A  l.p  ]  lbs 

1  s  dMg 


where  1^  is  steady  state  draft. 


10)  Plenum  Gage  Pressure 


(2, A. 7) 


(2. A. 8) 


(2. A. 9) 


pb  =  <pb  -  V 


lbs 

ft2” 


(2. A. 10) 


11)  Draft 


1*  =  Z  +  zc 
d  s 


ft 


(2. A. 11) 


The  system  parameters  and  constants  are  listed  as  follows 


1)  Adiabatic  Process  Coefficient 


y  =  1.4 
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2)  Leakage  Area 


A1  =  ZL  =  '438 


ft 


2 


3)  Leakage  Orifice  Coefficient 


C  =  .90 
n 


4)  Air  Density 


p  =  .002378 
Ka 


5)  Atmospheric  Pressure 


P  =  2116 

cl 


6)  Plenum  Area 


A,  =  200 

D 


ft 


2 


7)  Empty  Plenum  Volume 


V  =  383  ft3 

n 


8)  Craft  Weight 


W  =  6007 


lbs 
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9)  Area  of  Each  Immersed  Sidewall 
A  =  75/4  =  18.75  ft2 

10)  Draft  Initial  Condition 

L^Co)  =  .5  ft 

11)  Density  of  Water 

p  =  1.99 

12)  Gravity  Acceleration  Constant 

G  =  32.174  ft/sec2 

13)  Vertical  Center  of  Gravity  Location  Above  Keel 
Line 

Z  =  2.5  ft 

s 

14)  Number  of  Fans 

n  =  5 

15)  Steady  State  Fan  Output 

Q^(o)  =  35  ft3/sec 
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PRESSURE  LIFT  ONLY  SCHEMATIC  DIAGRAM 


B.  LINEAR  MODEL 


Prior  to  defining  the  system  states,  the  non  linear 
equations  were  linearized  around  their  steady  state  operating 
conditions  by  application  of  a  Taylor  series  expansion.  The 
details  of  the  linearization  process  is  shown  in  Appendix  A. 

The  resulting  simplified  pressure  only  signal  flow  graph 
for  the  linearized  system  is  shown  in  Figure  2. 

1 .  Transfer  Function  Development 

The  system  graphical  determinant  is; 

A  =  1  -  CL  +  L2)  +  L1  L2  (2.B.1) 

A  =  1  -  (L1  +  L2)  (2.B.2) 

where 

Lx  =  s'1s"1(-k8)(-k2)(k3)(-k4) 

=  -s‘2k2k3k4k8  (2.B.3) 

L„  =  s_1(k-)  (kJ(-kc) 

t  13  0 

=  -s^kgkgk^  (2.B.4) 


The  transmission  gain  from  the  control  input  to  the  accel¬ 
eration  output  for  calm  water  operation  is; 


LINEAR  PRESSURE 


Figure 


(2.B.5) 


K  = 


P1A1 


u 


where 


=  -k3k4k7k1Qs 


**  1 


(2.B.6) 


and 


A 


1 


1 


(2.B.7) 


and 


3 


w 


0 


C2.B.8) 


Substituting  equations  (2.B.2),  (2.B.6),  (2.B.7)  into 
equation  (2.B.5)  yields; 


-k3k4k7k10S 


-1 


u 


9"  =  - 5 - -| 

z  l-[-s'zk0k-,k„k0)  +  (-s"xk,krk7)] 


'2348' 


'3  6  7' 


(2.B.9) 


Reducing  further; 


_  "k3k4k7k10 

'9-  .  -j - 


u 


s  +k3kgk7  s+kjkgk^kg 


(2.B.10) 


The  transmission  gain  from  the  sinusoidal  sea  generator  input 
to  the  acceleration  output  with  no  control  signal  is; 


3 


fi 

z 


P2A2 


(2.B.11) 
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where : 


~klk2k3k4 


(2,B.12) 


and 


A 


1 


1 


(2.B.13) 


and 


3U  =  C  (2.B.14) 

Appropriate  substitution  into  equation  (2.B.11) 

yields 


-k1k2k3k4  3W 

l-[ (s"2k2k3k4kg )+ (-s-1k3kgk7 ) ] 


Reducing  further,  the  transfer  function  is 


3Z 

aw 


-klk2k3k4  s' 


s  +k3kgk7s+k2k3k4kg 


(2.B.15) 


For  the  condition  that  control  input  is  applied,  the 

acceleration  output  has  the  following  transform: 

-k.k^k-k..  _ s  3 
_  3  4  7  10  u 

g  -  y  1  r- ' 

z  s  +k3kgk7s+k2k3k4kg 

k  k  k  k  s2  3W 

-  -  (2.B.16) 

s  +k3kgk7s+k2k3k4kg 
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Following  the  same  procedure  for  the  applied  control 
input,  the  volumetric  plenum  air  mass  output  has  the 
following  transfer  function; 


A3  =  1  -  Lx 

(  2  ,  B  ,  1 7  ) 

P3  =  k10S 

(2.B.18) 

a2  =  1 

(2.B.19) 

P4  “  “klk2k3k6 

(2.B.20) 

y,  .  k10S+k2k3k4k8k10S  0 

m  s2+k3kgk7s+k2k3k4k8  u 

(2 ,B.21) 

Equation  (2,B,15)  is  of  the  form; 

OUTPUT  -as2 

INPOT"  “  „2,x  . 

S  +DS+C 

( 2 .B. 22) 

and  with  no  control  action,  becomes; 

y(s)  _  as2 

(2 .B. 23) 

where 

yCs)  =  3Z 

(  2  ,  B  .  2  4  ) 

and 

2V(s)  V 

WVsT  - 

C2.B.25) 
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where  the  coefficients  consist  of  the  following  parameters; 


a  -= 


-k1k2k3k4 


b  = 


k3k6k7 


k2k3k4k8 


d  => 


~klk2k3k6 


(  2  ,  B  ,  2  6  ) 

(2.B.27) 

(2.B. 28) 

(2 .B. 29) 


The  values  of  the  coefficients,  k  ,  for  the  XR-3  craft  are 
listed  in  Table  I  and  were  obtained  from  the  results  of  the 
linearization  process  of  Appendix  A. 

Special  attention  should  be  given  to  the  coefficient 
that  contains  k^  which  is  a  sine  function  of  the  form: 

2W  V  w .  L 

k  =  [  (  — E__W)  sin  (*£-£) ]  (2.B.30) 

1  a).  2Vw 


k^  may  be  simplified  by  an  approximation  detailed  in  Appendix 
A  of  Ref.  8.  The  simplification  is  based  upon  the  fact  that 
for  ahead  seas,  the  incident  frequency,  uk  is; 


Cl) 


i 


1- 


(2.B.31) 


For  convenience  in  this  study,  the  value  of  the  ships  velocity, 

Vg ,  is  set  equal  to  the  gravitational  constant,  G.  This  value 

of  V  ,  when  converted  to  knots,  is  19.3  knots  which  is  within 
s  ’  ’ 
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TABLE  I 

SYSTEM  COEFFICIENT  (k  )  VALUES 


n  FORMULA 

2WpVw  '  W,  Lp 

i  W 


V57 


XR~3  VALUE 
2267 ,782 


"7^2T8'3'lL9y  • 


Jw 


3 ,5336x10 


~3 


:3  yPb(O) 


2996.0 


Ab 

M 


1.0711 


nk 


C  A  / 

+  -4-i  / 


2Py(0) 


6.14422 


'7  V  (0) 
m 


3 .5052x10 


-3 


k8  Ab 


200 . 0 


k1Q  CONSTANT 


1.0 
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the  operating  range  of  the  XR-3,  Continuing; 

1-  /l +4U 


w , 


c: 


-2 


It  follows  that 


0) .  =0) 
l  e 


(J 


This  approximation  is  valid  for  all  coe>2,5  radians/sec. 
the  definition  of  the  incident  wave, 


2  _  2TTg 


C 


k^  becomes 


«iLp 


Sin'  2 V 


Cd 


0.  L 


2  V 


w 


c 


Where  the  horizontal  area  of  the  plenum  chamber  is 


Afa  =  WpLp 


C 


For  a  speed  of  19.3  knots  and  since 


JW 


2j£ 

w 

e 


c: 


,B. 32) 


.B. 33) 

From 


.B. 34) 


.  B  .  3  5  ) 


.  B  ,  36  ) 


.  B  .  37  ) 
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kn  simplifies  to j 


For  the  XR-3  , 

A,  =  200  ft2 
b 

L  =  20  ft 
P 


Then 


a  = 


2267,782 


e.  ,62.8319. 
Sm('  ■  v'"’  • »  > 

LW 


62 .8319, 


From  Equation  (2.B.38),  it  can  be  seen  that  a  null 
will  occur  for  nL^  =  and  a  peak  value  for  k^ 
for  nL^  =  L^/2  where  n  is  an  integer. 

2 ,  State  Equation  Development 

The  state  equations  are  of  the  form 


X  =  AX  +  B3u  +  G3W 


where  the  state  variable  vector  is  chosen  as 

X  =  (3  ,  3*  ,  3 V  )T 
—  z  ’  z  ’  m 


and  the  output  variable 


y  =  3  =  C  X  +  F  3W 

Z  N  — 


(2.B.38) 


(2.B.39) 

value  for 
will  occur 


(2.B.40) 

(2.B.41) 


(2.B.42) 
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The  states  are  defined  as 


XI  =  3Z  =  (Z-Z  ),  Incremental  Heave(Draft)  C2.B.43) 
ss  ' 

X2  =  3Z  =  (Z~Z  ),  Incremental  Heave  Rate  C2.B.44) 

S3  ’ 


X3  =  3V  s  (V  "V  ) j  Incremental  Air  Mass 
m  m  mss  '  -- 


Volume 


where  V  =  t— 
m  p  „ 


(2.B.45) 


(2.B.46) 


with  control  defined  as  a  scaler,  rate  of  incremental  change 
of  air  mass  volume. 


3  =  3V 

u  m. 


(2.B.47) 


and  where  the  sea  input  disturbance  is  3W. 

From  the  signal  flow  graph  of  Fig.  2, 


3V  =  (3  k  =kc3P,)s 
m  u  10  6  b 


-1 


(2.B.48) 


X1  =  32  *  5-1  *= 


(2.B.49) 


9  •  —1  «, 

X,  =  3  =  s  x  3z 

1  z 


C2.B.50) 


X2  =  3z  =  "  3Pbk4 


( 2 , B , 51 ) 
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where 


3P  =OV  k_-k_[-3Wk,-kQ  3Z])k, 
b  m  7  2  18  3 

Rearranging 

3P,  =  3V  k-k-H-k-k-k- 3W+k„k0k03Z 
t)  m  O  /  J.  2  0  zoo 

and  grouping  in  order  of  states 


C2.B.52) 


(2.B.53) 


3P,  =  k.k.ko3Z+k0k„3V  +k^k^k^3W 
D  2  o  o  6  I  jn  12  3 

The  output  term,  y'is 


(2.B.54) 


y  = 


(2.B.55) 


Substituting  equation  (2,B,5)  into  equation  (2.B.55)  yields, 


y=-k-k_k.  kQ  3Z-k_k..k-  3V  -k,  kok0k.,  3W 


(2.B.56) 


Finally,  the  incremental  change  in  the  volumetric  air  mass 
is  written  as 


X3=3Vm=<5uk10-k6 9Pb> 

“  -k2k3k6k83Z-k3k61'79Vm 
k^kjhgkg  3W+k^Q  3u 


( 2 , B , 57  ) 


(2.B.58) 
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Equations  (2.B.50),  (2,B,56)  and  (2,B.57)  are  the  states  of 
interest  and  are  grouped  below, 


XI  =  X2 


X2  =-k2k3k  k8xi«k3k4k7-X3-.k1k2k3k43W 


y  =  -k2k3k4kgXl-'k2k4k7X3-k^k2k3k4  <W 


Thus  the  linear  system  matrices  are: 


A  =  -k2k3k4kg  0  ->kgk4ky 

_k2k3k6k8  0  "k3k6K7 


(2.B.64)  F=-k1k2k3k4 


G  =  -k1k2k3k4 


-k^kjkgkg 


(2.B.59) 


(2.B.60) 


X3  =  -k2k3k6kgXl-k3k6k7X3-k1k2k3k6 3W+k103u  (2.B.61) 


with  the  output  acceleration: 


(2.B.62) 


(2 -B.  63) 


(2.B.65) 


(  2  ,  B  .  6  6 ) 
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Computer  simulation  of  the  linear  system  was  accom¬ 
plished  using  the  Naval  Postgraduate  School  program  "IODE  : 
operated  from  the  CPCMS  terminal  of  the  IBM-360.  See  Appendix 
C  for  IODE  listing.  The  program  ocean  generator  for  3W  inputs 
a  single  frequency  at  a  time  using  the  special  function  fea¬ 
ture,  Wavelength  (Ly)  will  change  with  encounter  frequency 
<w  )  as  shown  in  equation  (2.B.35),  It  should  be  remembered 
that  the  craft  speed  is  constant  at  19,3  knots  for  the  data 
base  of  this  design, 

3 .  The  Characteristic  Equation 

The  characteristic  equation  under  no  control  condi¬ 
tions  of  operation,  is  formed  from  the  A  matrix  (see  Appendix 
A), 


DET(SI-A) 


(2.B.67) 


DET 


s  -10 

2267.782  s  11.248 

13002.8  0  64.493s 


(2.B.68) 


Then  the  characteristic  equation  is; 

s3+64.493s2+2267,782s+,5701  (2.B.69) 

which  factors  into 

(s  +  2, 5414xl0’'4)(9+3  2. 2  5+J3  5.0  4KS  +  3  2. 2  5-j  3  5. 04  ) 

C2.B.70) 
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and  is  shown  in  the  pole^zero  map  of  Fig.  3,  Note  that  the 
dominant  pole  is  a  real  pole  close  to  the  origin.  Also  note 
that  the  complex  pair  will  affect  early  time  response  since 
they  are  relatively  far  removed  from  the  origin, 

A  root  locus  plot  with  no  control  input  is  ?hown 
in  Fig.  4.  The  plot  covers  a  range  of 

5  <  kg  <  26  (2.B.71) 

Table  II  lists  the  characteristic  equation  roots  for  the 
range 


5  <  kg  <50 


(2.B.72) 


As  can  be  seen  from  equation  (2.B.63),  kg  is  present 
in  both  the  a31  and  a33  system  coefficients.  In  order  to 
perform  the  root  locus  study,  the  characteristic  equation 
was  rearranged  to  show  the  adjustable  parameter  kg ; 

3  2  2 

s  +k3k7s  +k25<3k4kgS+kg  (k3k4k7k2k3kg+k2k3k4kgk3k7  )  +  kgS 


( 2 . B . 73a) 


with  parameter  values 


s3+(kg+10, 5016) s2+2267 ,782 s+ ,0928kg 


( 2 , B . 73b) 


From  the  table  and  root  locus  plot,  it  can  be  seen  that  as  kg 
becomes  large,  the  complex  root  pair  become  real  with  one 


root  moving  toward  the  origin  and  the  other  migrating  to 

negative  infinity.  Additionally,  the  pole  near  the  origin 

moves  very  little  and  remains  real  throughout  the  excursion 

of  kg.  It  was  observed  that  this  was  very  similar  to  the 

root  locus  plot  in  the  air  flow  study  conducted  in  Ref,  1. 

However,  the  removal  of  the  buoyancy  term  in  this  study 

caused  the  pole  near  the  origin  to  stay  real  whereas  the 

same  pole  in  Ref.  1  became  immaginary  for  large  values  of  k_ . 

b 

The  root  locus  demonstrated  the  high  sensitivity  of  the  com¬ 
plex  pair  to  changes  in  air  flow  rate ,  while  the  pole  near 
the  origin  exhibited  a  very  low  sensitivity  to  changes  in 
kg.  This  observation  further  demonstrated  that  the  vertical 
oscillatory  characteristics  of  the  ship  depend  upon  changes 
in  kg . 

As  a  further  check  on  the  effect  of  air  flow  on 

C.G.  acceleration,  an  investigation  of  the  effect  of  changing 

kg  on  the  frequency  spectrum  of  C.G.  accleration  was  conducted 

and  the  results  are  shown  in  Appendix  D.  From  these  results 

it  seemed  appropriate  to  seek  a  control  law  that  would  adjust 

ke ,  the  air  flow  rate.  In  an  effort  to  gain  a  better  under- 
b 

standing  of  the  effect  of  various  air  flow  rates  on  the  sys¬ 
tem,  a  sensitivity  and  stability  analysis  was  performed  in 
section  IVB,  The  ROOTLO  program  used  in  the  above  study  is 
shown  in  Appendix  C, 
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X-SCfllE=2. 00E+0 1  UNITS  INCH. 
Y-SCflLE=2.00E+01  UNITS  INCH-. 


ROOT  LOCUS  PLOT  FOR  MODIFIED 
PRESSURE  ONLY  LINERR  MODEL 


Figure  4 
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Ill,  VERIFICATION  OF  THE  MODEL 


A.  LINEAR  AND  NONLINEAR  SYSTEM  VALIDATION 

The  validation  of  the  simplified  model  of  a  single  vari- 
ble  pressure  lift  force  was  conducted  by  comparison  of  its 
frequency  response  characteristics  with  those  of  the  model 
presented  in  Ref.  1,  Several  techniques  were  used  to 
accomplish  this, 

1.  Single  frequency  constant  amplitude  sinusoids  were 
applied  to  both  the  nonlinear  system  of  Fig.  1  and 
also  the  linear  system  of  Fig.  2.  The  acceleration 
output  signal  for  the  nonlinear  system  is  shown  in 
Fig.  5.  Its  sinusoidal  form,  in  steady  state  indi¬ 
cates  the  linearity  of  the  system,  which  was  the  case 
throughout  the  frequency  range. 

2.  The  linear  transfer  function  of  center  of  gravity 
acceleration  (Eq.  2.B.15)  was  subjected  to  a  constant 
wave  height  input  using  a  spectrum  of  discrete  en¬ 
counter  frequencies.  The  heave  accelerations  were 
converted  to  dB  and  plotted  in  a  bode  diagram  shown 
in  Fig.  6 

3.  Discrete  fourier  transform  technique,  developed  in 
an  independent  study  by  S,  Carpenter,  [Ref,  9],  was 
also  used  to  test  the  response  of  the  linear  system 
(Fig.  9)  to  an  irregular  sea  input  shown  in  Fig.  7 
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consisting  of  twelve  components  of  frequency  and 
amplitude  from  the  Pierson-Moskowitz  energy  density 
spectrum  given  in  Fig.  8  [Ref.  8,  9], 

The  results  of  procedures  2  and  3  above  were  compared 
with  the  linear  system  response  of  Ref.  1,  The  magnitudes 
for  the  center  of  gravity  acceleration  showed  approximately 
20%  difference  between  the  pressure  lift  only,  no  control, 
response  curve  (Fig,  6)  and  the  linear  system  of  Ref.  1 
(Fig.  10).  The  general  shape  of  the  curves  are  identical. 
Therefore,  the  pressure  lift  only  model  was  assumed  to  be 
valid  and  useful  for  linear  regulator  control  design.  The 
observed  peaks  and  nulls  of  Figs.  6,  10  are  a  result  of  the 
sine  function  numerator  coefficient,  k^ ,  as  described  in 
Section  III.  The  relationship  between  plenum  length  (L^) 
and  wavelength  (L^)  expressed  in  k^ ,  determine  where,  in 
the  encounter  frequency  spectrum,  the  nulls  and  peaks  occur. 
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Figure  7 
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Figure  8 
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Figure  9 
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ENCOUNTER 


IV,  THE  CONTROL,  LAW 


A.  LINEAR  REGULATOR  DESIGN  APPROACH 

The  linear  pressure  lift  only  model  is  described  in 


state  equation  form; 

X1 

0  10 

X1 

0 

X2 

= 

a21  0  a23 

X2 

+ 

0 

*3 

a31  0  a3  3 

X3 

1 

— .  — 

_  _ 

—  — 

with  the  cost  function  of  the  form 


J 


(XTQX  +  auTB3u) 


(4. A.  2) 


where  the  final  value  error  is  not  included.  The  system  may 
be  designed  using  the  linear  regulator  technique  with  the 
given  quadratic  cost  function.  An  optimal  control  was 
found  using  both  the  transient  and  steady  state  solutions 
of  the  Riccati  equation.  The  control  law 


*  -IT 
3U  =tR  B  KX 


(4. A. 3) 


where  the  K  matrix  was  found  from  the  solution  of  the  Riccati 
equation 
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(4, A. 4) 


K  =  -KA  -  ATK  -  Q  +  K  B  K''1BTK 

-  rs  -<  '  1  J 


with  the  appropriate  choice  of  the  Q  and  R  coefficients,  the 
feedback  gains  for  the  system  were  calculated.  Unfortunately, 
the  first  choice  for  Q  and  R  may  not  yield  control  gains  which 

*■*  W 

are  physically  realizable  or  practical.  For  this  reason,  it 
is  appropriate  to  examine  the  effects  of  the  control  gains  on 
the  pole  location  in  the  s  plane  for  the  cor- trolled  system. 
Such  a  study  was  done  and  is  presented  in  part  B  of  this 
section.  The  choices  for  Q  and  R  were  derived  by  trial  and 
error.  However,  the  design  criterion  was  to  reduce  this  mag¬ 
nitude  of  the  frequency  response  curves  through  the  study  of 

various  Bode  plots.  To  this  end,  many  Q  and  R  values  were 

?*<  ^ 

tried . 

The  Riccati  equation  computations  were  done  using  a  com¬ 
puter  program  documented  in  Ref,  10  and  modified  by  B. 
DesJardins  in  Ref.  11,  Appendix  C  illustrates  the  input/ 
output  format  of  this  program.  Based  upon  the  above  design 
criterion,  the  final  weighting  matrix  chosen  was; 


Q  = 


0 


0 


0 


and 


R  =  .03 


(4. A, 5) 
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The  optimal  control  gains  obtained  using  the  above  Q  and  R 

rs 

matrixes  were ; 

gxl  =  -1 , 689xl03 

g12  =  1 , 19xlQ2  (4. A. 6) 

g13  =  -1.83X101 

These  gains  were  multiplied  by  their  respective  states  and 
applied  to  the  state  equations  (2.B.59,  61,  62).  Under  state 
feedback  control,  the  resultant  output  acceleration  frequency 
magnitude  was  reduced  by  13%  from  the  peak  no  control  magni¬ 
tude.  Table  III  lists  the  peak  lobe  values  of  the  C.G. 
Acceleration/wave  height  magnitudes  versus  frequency  for  the 
nonlinear  and  linear  systems.  Each  system  was  exercised 
under  control  and  no  control  conditions.  The  steady  state 
sinusoidal  output  for  the  linear  model,  under  control  and 

no  control  conditions  is  shown  in  Table  IV  for  u>  =45  rad/sec, 

e 

where 


3u  =  gllXl  +  S12X2  +  S13X3  (4. A. 7) 

and 

y  =  9Z  C4.A.8) 

The  control  signal  for  the  nonlinear  system  was  generated  by 

using  the  linear  control  low  coefficient  found  using  linear 

regulator  theory.  The  reduction  in  acceleration  was  not 

uniform  throughout  the  frequency  range.  At  law  frequencies, 

we<l’7  rad/sec,  control  gains  caused  an  increase  in  C.G.  acceleration. 
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As  uie  approached  100  rad/sec ,  C,G,  acceleration  approached  but 
never  exceeded  the  no  control  condition.  Effective  control 
was  realized  in  the  frequency  range  17<We<100  radians/ sec. 


B.  SENSITIVITY  .AND  STABILITY 

Since  a  set  of  optimal  gains  were  found  for  the  system, 
it  remained  to  test  the  system  under  feedback  control  and 
observe  the  sensitivity  and  stability  for  acceptability. 

The  locus  of  roots  of  the  characteristic  equation  for 
the  pressure  lift  only  system  without  state  feedback  control 
was  shown  in  Section  III,  The  objective  of  state  feedback 
control  is  to  adjust  the  position  of  the  complex  pair  of 
roots  in  such  a  way  as  to  reduce  center  of  gravity  accelera¬ 
tions.  For  the  open  control  loop,  the  transfer  function 
was  shown  (Eq.  2,B,22)  to  be 


•  • 

3Z  . 

W  ~ 


-as 


s  +bs+c 


(4.B.1) 


Addition  of  state  feedback  will  change  the  value  of  b  and  c. 
Note  that  if  the  k^  sj.nc  function  null  effects  are  removed 
from  the  numerator ,  then  a=c .  It  follows  that  the  magnitude 
of  the  Bode  plot  will  reach  unity  value  at  high  frequency 
regardless  of  changes  in  b  and  c,  In  other  words,  the  sys¬ 
tem  has  high  pass  frequency  characteristics.  Thus,  the 
results  of  increasing  the  value  of  c  will  be  a  reduction  of 
the  magnitude  plot  in  the  lower  range  of  frequencies  of 
interest.  Also,  note  that  changes  in  the  b  coefficient  will 
affect  the  resonant  peaking  at  the  corner  frequency. 


TABLE  III 

PEAK  LOBE  VALUES  OF  FREQUENCY  RESPONSE,  PRESSURE  LIST  ONLY 
LINEAR  S  NONLINEAR  MODELS  UNDER  CONTROL  AND  NO  CONTROL 

CONDITIONS 


u  NONLINEAR  LINEAR 

e 


NO  CONTROL 

CONTROL 

NO  CONTROL 

CONTROL 

6 

18,6 

26,5 

18,5 

27,8 

7 

18,5 

25,2 

18.6 

26,6 

8 

15.7 

20,3 

15,7 

21.4 

9 

9.7 

12 . 7 

9.8 

12 . 7 

14 

42.2 

42.2 

42.3 

45.7 

15 

48.2 

48 . 2 

48.4 

50.9 

16 

49.7 

49 . 7 

49,9 

51.1 

17 

45.9 

45,9 

46.1 

46.2 

24 

69.8 

62.9 

70 . 5 

63.4 

25 

78.1 

69,7 

79.0 

70.3 

26 

79.1 

70.0 

79.9 

70.6 

27 

72.1 

63,4 

72.8 

63.9 

28 

34 

90.8 

77.7 

91.8 

78.5 

35 

100 . 8 

86.3 

101.8 

87 . 3 

36 

102.4 

86 . 8 

102.3 

87.9 

37 

91.7 

78.6 

92.6 

79.6 

38 

44 

99.5 

86.8 

100.1 

87,4 

45 

110.0 

96.4 

113.0 

98.6 

46 

110,0 

96,8 

111.5 

98 . 0 

47 

99.7 

88,0 

101.0 

88.9 

48 

54 

97.1 

88.5 

97,5 

88.5 

55 

107,7 

98,5 

108.0 

98,7 

56 

108.1 

99,2 

108.8 

100,0 

57 

98.2 

90.5 

99.0 

91,2 

TABLE  XIX 
(continued) 


co  NONLINEAR  LINEAR 

e 

NO  CONTROL  CONTROL  NO  CONTROL  CONTROL 


64 

88.7 

83,8 

89.3 

84.3 

65 

99.3 

94,1 

99.5 

94.8 

66 

100.4 

95,4 

101.1 

96,5 

67 

91,9 

87.6 

91.9 

88,0 

68 

69 

74 

79,0 

76.5 

79,4 

77,0 

75 

89.5 

86,9 

89.9 

87 . 3 

76 

91.3 

88 . 8 

91.3 

88.6 

77 

84.4 

82,2 

84.8 

82.6 

78 

79 

84 

69,9 

68.8 

70.1 

69,0 

85 

80.2 

78.9 

80.1 

79.3 

86 

82.6 

81.4 

83.0 

81.7 

87 

77,4 

76.1 

77.0 

76.5 

94 

61.9 

61.4 

62.0 

61.5 

95 

71.9 

71,5 

72.1 

71.6 

96 

75.0 

74 . 5 

75.2 

74.7 

97 

70.9 

70.5 

70.9 

70.6 
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TABLE  IV 


SINUSOIDAL 

LINEAR  MODEL 
AT 

OUTPUT  WITH/WITHOUT 

to  -45,0  Rad/Sec . 
e 

FEEDBACK  CONTROL 

TIME 

CONTROL 

CONTROL  EFFORT 

NO  CONTROL 

T 

Y 

9 

u 

Y 

0.780 

54,38 

4,28 

61.02 

0.785 

34.33 

2 . 64 

38 . 03 

0.790 

12,50 

0.86 

13 . 08 

0.795 

-9.96 

-0,95 

-12 . 54 

0.800 

-31,91 

-2,72 

-37.52 

0.805 

-52.23 

-4,35 

-60.56 

0.810 

-69,85 

-5 .76 

-80.49 

0.815 

-83,86 

-6.87 

-96.26 

0.820 

-93.56 

-7,62 

-107 . 08 

0.825 

-98.44 

-7,99 

-112 .38 

0.83  0 

-98.23 

-7.94 

-111.88 

0.835 

-92.97 

-7.48 

-105.62 

0.840 

-82 . 90 

-6.64 

-93.92 

0.845 

-68 .60 

-5,45 

-77.37 

0.850 

-50.74 

-3 .97 

-56.85 

0.855 

-30,27 

-2.32 

-33.39 

0.860 

-8 . 34 

-0,53 

-8 . 21 

0.865 

14.22 

1.29 

17.40 

0.870 

35.95 

3.04 

42.10 

0.875 

55.82 

4.63 

64.64 

0,880 

72.82 

5.99 

83 .85 

0,885 

86,06 

7.03 

98.73 

0.890 

94.87 

7.72 

108,53 

0.895 

98 . 79 

8 . 00 

112.73 

0.900 

97,61 

7.88 

111.11 

0.905 

91,42 

7  .  34 

103.79 

0.910 

80,50 

6.43 

91,11 

0.915 

65.44 

5.19 

73.74 

0.920 

47.01 

3,68 

52.56 

0.925 

26,49 

1.98 

28.68 

0.930 

3,94 

0.17 

3  .  31 
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The  Characteristic  Equation  For  Feedback  Control 


The  signal  flow  graph  of  Fig.  2  was  modified  to  show 
feedback  gains  and  is  shown  in  Fig,  11.  Using  Masons 
gain  rule,  the  characteristic  equation  was  derived  from  the 
S.F.G.  (Eq.  4.B.10),  The  resulting  equation  was  compared 
to  characteristic  equation  derived  from  the  modified  system 
matrix  with  feedback  of  Appendix  D.  It  was  noted  that  they 
differ  by  the  terms 


a21a3  3  and  ->a2  3^3  • 

T;hese  terms  very  nearly  cancel  and  are  insignificant  (a^ 
a33-a23a3i= • 570D  when  compared  with  the  other  larger  terms 
in  the  characteristic  equation.  Therefore,  their  omission 
from  the  signal  flow  graph  of  Fig.  11  is  valid.  The  deri¬ 
vation  follows : 


~§3  _  “S3 


1  s 


(4.B.2) 


(kg )  ( -kg ) ( -ky ) 


"k3k6k7 


(4.B.3) 


(-ky  )  (kg)  (-k^)  (-kg)  '"kykgk^kg 


(4.B.4) 


(k?)(k3)(-k4)Cg2) 


'k3k4k7 


g2  ( 4 , B . 5 ) 


(ky)(k3) (-k4) (-g1)  k3k4ky 


g  ^  (4  ,B,  6) 


i 


Figure  11 

Linear  Pressure  Lift  Qnly 
Signal  Plow  Crapli  With  Stale  Feedback  Control 


4  ;  l-[L1tL2tL3  VV  ♦  Ll  L3 


(4.B.7) 


A  =  l^C- 


l33  a21  a23  .  a23 

„  r  g2  +  _T 

s  s  s 


g^K^-X. 


-a 


21- 


(4.B.8) 


A  = 


(s3  + 


S3S 


*  a33s2  * 


a21S  +  a23 


g9s 


a23Sl 


+ 


a21S3  } 


(4.B.9) 


s3+(a33+g: 


)s*  +  (a 


23' 


:2+a21)s 


+  (a 


21s3 


23sl 


) 


(4.B.10) 

2 .  Root  Locus 

The  following  root  locus  plots  are  a  study  of  the 

effects  of  changing  the  three  different  control  feedback  gains. 

Many  weighting  combinations  were  tried  and  it  was  discovered 

that  the  system  was  most  sensitive  to  changes  in  air  flow  rates. 

This  observation  was  in  agreement  with  the  analysis  of  Section 

II  that  the  system  was  most  sensitive  to  k. .  Each  of  the 

b 

following  root  locus  plots  and  associated  table  of  values  was 
chosen  as  representative  of  the  alternatives  tried.  The 
characteristic  equation  used  is  shown  as  Eq,  (4.B.10).  The 
Naval  Postgraduate  School  program  for  the  IBM  360/67  computer 
was  used  to  generate  the  tables  and  plots.  The  program 
listing  is  shown  in  Appendix  C, 
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In  the  first  study,  the  constant  gains  used  were 


S11  =  425,0  »  §i2  =  91 

and  g13  was  varied  over  three  decades  starting  from  g^3  =  .1. 
The  resultant  plot  and  table  of  values  are  shown  in  Fig.  12 
and  Table  IV  respectively.  For  ease  of  reference,  the  sys* 
tem  poles  were  given  letter  designations  for  the  characteris¬ 
tic  equation 

(S+A) (S+B) (S+C)  =  0  (4.B.11) 

where 

A  =  (+32.25-j35.04) 

B  =  (+32.25+j35.04) 

C  =  (  2 . 5xl0""4  + j  0  ) 

As  g^3  nears  -20,  pole  a  rapidly  enters  the  real  axis  and 
migrates  towards  negative  infinity  and  pole  B  moves  upward 
across  the  real  axis,  Pole  C  moves  in  a  negative  direction 
until  g13=-26  and  then  moves  below  the.  real  axis  tracking  the 
mirror  image  path  of  pole  B.  It  was  seen  from  this  plot 
that  as  g13  nears  -20  that  the  complex  pair  undergoes  the 
greatest  movement.  In  this  study  there  was  no  point  at 
which  instability  occurred,  When  the  gains  were  tested  on 
the  system,  there  was  no  noticable  decrease  in  vertical 
accelerat ion . 

The  next  phase  of  this  investigation  uses  constant 

gains  of 
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X-SCfll_E  =  2.  OOE  +  0 1  UNITS  INCH. 
Y-SCRLE=2. OOE+Ol  UNITS  INCH. 


ROOT  LOCUS  PLOT  FOR  MODIFIED 
PRESSURE  ONLY  LINEAR  MODEL 


Figure  12 


TABLE  IV 

Roots  For  Modified  Pressure  Lift  Only  Model, 

gl3  Varies 


^6 

A 

B 

-0,099 

-31 . 49+ j46 . 96 

-31.49-546 , 96 

-0,127 

-31,52+j46.88 

—31 . 5  2  — j  4  6 . 8  8 

-0.162 

-3 1 . 5  2  + j  4  6 , 8  6 

-31 . 52-j  46.86 

-0.207 

-31.53+j46.84 

-31,53-j46,84 

-0.263 

-31 , 53  + j  46 , 81 

-31 . 5  3— j  4  6 . 81 

-0,335 

-31.54+546,76 

-31.54-j46.76 

-0.427 

-31 . 55  +  j  46,71 

-31. 55-j 46,71 

-0,544 

-31 , 57  + j  46 , 64 

— 31,57— j46, 64 

-0.693 

-31,59+j46,56 

-31.59-j46 . 56 

-0 ,883 

-31 . 61+ j  46.45 

-31.61-j46,45 

-1.125 

-31 . 64+ j  46.30 

-31 , 64- j  46 , 3  0 

-1.433 

-31,68+j46.12 

-31. 6 8-j  46 .12 

-1.825 

-31.7  2_j  45,88 

-31.72-j45.88 

-2.326 

-31 . 77  + j  45 . 58 

-31,77-345.58 

-2.962 

-31 , 84+ j  45.18 

-31 . 84-j  45.18 

-3.773 

-31. 91_j44 .67 

-31 , 91-3  44.67 

-4 .807 

-31, 99  +  j  4  3 .99 

-31.99-j43.99 

-6.124 

-32.06+j43,08 

-32 . 06- j  43 , 08 

-7.801 

-32. 09+ j 41, 87 

-32 , 09-j41, 87 

-9.937 

-32 . 03  +  j  40 . 20 

-32.03-j40.20 

-12.658 

-31 . 70  + j  3  7 . 8  7 

-31.70-j37,87 

-16.125 

-30. 60+j34, 50 

-30.60-j34.50 

-20.540 

-26.95+j30,22 

-26.95-j30.22 

-26,165 

-48 . 97+j  0 

-20 . 53  + j  29 . 31 

-33.331 

-64. 98+jO 

-16.02+j30,86 

-42.458 

—  80.2  6  + j  0 

-12.84+j32.65 

-54.085 

-96 , 5 5+ j  0 

-10 , 36  +  j  34 , 38 

-68.90 

-115 , 02  +  j  0 

-8 ,36+36.02 

^r425-0-  s12 

C 

-1.49+jO 
-1 . 586  + j  0 
-1,611+jO 
-1 . 64 3  + j  0 
-1 , 684  + j  0 
-1. 736+j 0 
-1, 802+j0 
-1.887+jO 
-1.995+jO 
-2 . 134  + j  0 
-2,312+jO 
-2 . 540  + j  0 
-2 , 835  +  j  0 
-3 . 216+ j  0 
-3.711+jO 

-4. 3 5 9  +  j  0 
-5 . 214  + j  0 
—  6.359  +  j0 
-7,923+jO 
-10 . 134  + j  0 
-13 . 450  + j  0 
-19 . 0  3  2  + j  0 
-30 . 645  + j  0 
-20 . 534-j  29, 31 
-16 . 0 23- j 30 . 86 
-12 , 83  8-j  32,65 
-10 . 367-j  34 , 38 
-8 . 3  6  3-j 36 . 02 
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"  "195.32,  g^2 


34 .98 


while  varied  over  two  decades  from  -.1.  These  gains 
were  calculated  from  the  Riccati  with  weighting  matrix; 


Q  = 


10  0  0 

0  122  0 

0  0  .5 


and 


C4.B.12) 


R  =  .03 


(4.B.13) 


The  resulting  plot  and  table  of  pole  values  are  shown  in 
Fig.  13  and  Table  V  respectively.  As  g-  varies,  the  poles 
follow  similar  paths  as  in  the  previous  study.  The  greatest 
degree  of  pole  movement,  or  sensitivity,  occurs  at  g^g=-16.0. 
This  set  of  gains  produced  instability  when  g^g>-1.08. 

In  the  final  study,  the  constant  gains  were 


S11  =  -1688 . 095  ,  g  2  =  119.16,  g±  =  -18.306 

Again,  g^g  was  made  to  vary  over  three  decades  starting  with 
g^g  =  -.1.  The  resultant  plot  and  table  of  values  are  shown 
in  Fig.  14  and  Table  VI  respectively.  These  gains  were  cal¬ 
culated  using  the  Q  and  R  values  shown  in  (4,A.5).  Again, 
the  poles  follow  similar  paths  but  pass  closer  to  the  origin 
than  in  the  previous  two  studies,  causing  a  greater  net 
decrease  in  the  systems  oscillatory  tendencies.  As  in  the 
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X-SCRIE=2.00E+Q1  UNITS  INCH. 
T -SCALE=2 . OOE  +C 1  UNITS  INCH. 


ROOT  LOCUS  PLOT  FOR  MODIFIED 
PRESSURE  ONLY  LINERR  MODEL 


Figure  13 


TABLE  V 


Roots  For 

Modified  Pressure 
g12=34r98 

Lift  Only  Model, 
,  g13=varies 

g.^-195 . 32  , 

^6 

A 

B 

C 

— 

— 

— 

-0.099 

-32 , 65+j40 ,60 

-32,56-j40,60 

0, 81  + j  0 

-0.127 

-32, 66  + j  40 , 50 

-32 , 66-j40 , 50 

0 . 71+ j  0 

-0.162 

-32.67+j40.48 

-32.67-j40,48 

0 . 6  8  +  j  0 

-0,207 

-32 . 67+j4G , 45 

-32,69-j40,45 

0 . 64+ j  0 

-0,263 

-32 . 67  + j  40 . 41 

-32.67-j40,41 

0 . 6  0  + j  0 

-0.335 

-32 . 68  +  j  40,35 

-32 . 68-j  40 , 3  5 

0 . 54  + j  0 

-0.427 

-32 , 69  +  j40 . 29 

-32 . 6 9- j  40 . 29 

0 . 46  + j  0 

-0.544 

-32,7  0+ j  40.20 

-32.70-j40.20 

0 . 37  + j  0 

-0.693 

-32,71+ j40, 09 

-3  2 . 71- j  40 . 10 

0 , 2  5  + j  0 

-0.883 

-32,72+j39.95 

—  32 .72  —  j39.95 

- . 09  + j  0 

-1.125 

-32.74+j39,78 

-32 , 74—j 39.78 

-Q.ll+jO 

-1.433 

-32,76+j39.54 

-32 . 76- j  3 9 . 54 

-0. 37+j0 

-1,826 

-32 . 78  + j  39.25 

-3  2. 7 8- j  3  9.24 

-0. 71+j0 

-2 .326 

-32.81+j38.86 

-32,81-j38.86 

-1.14+j0 

-2,962 

-32 , 84  + j  38.35 

-32 . 8  4 — j 38.35 

-1.71+j0 

-3.774 

-32.86+j37.68 

-32 . 8  6  —  j 37.68 

-2.46+jO 

-4.807 

-32 . 86  + j  36,79 

-32.86-j36.79 

-3 . 47  + j  0 

-6.124 

-32 , 81  + j  35,59 

-32 , 81- j  35,59 

-4.85+jO 

-7.801 

-32.66+j33. 93 

—  32 . 66— j33 . 93 

-6 . 8 0  + j  0 

-9,937 

-32.23+j31.54 

-32.23-j31.54 

-9 . 74  + j  0 

-12.66 

-31 . 01+ j  27.93 

-31 . 01- j  27.93 

-11.48+jO 

-16.12 

-26.18+j22. 73 

-26.18-j22.73 

-27 , 87  + j  0 

-20.54 

-47 , 81+ j  0 

-18 . 37-j  23.83 

-18. 37-j 23. 83 

-26.16 

-61 . 49  +  j  0 

-14 , 27-j  26 . 51 

-14 , 27-j  26 . 51 

-33,33 

-74 . 16  +  j  0 

-11 , 43-j  28 . 89 

-11 . 43-j  28 . 88 

-42,46 

-8  7 , 51+j  0 

-92,12-j31, 05 

-9 , 21- j  31.05 

-54.09 

-102 , 51+ j  0 

-7 3, 88 -j 33. 04 

-7 , 39- j  3  3 , 04 

-68.90 

-120. OO+jO 

-58 . 68-j  34 , 90 

-5. 87- j  34 , 90 

X-SCRLE=2 . OOE+O 1  UNITS  INCH. 
T-SCRLE=2.00E+01  UNITS  INCH. 


ROOT  LOCUS  PLOT  FOR  MODIFIED 
PRESSURE  ONLY  LINEAR  MODEL 


Figure  14 


TABLE  VI 

Roots  For  Modified  Pressure  Lift  Only  Model,  g  =-1688,  g  .=119, 

gl3=^18 , 3(varies)  1 


^6 

A 

B 

C 

-0.099 

-34 , 66  + j  52,35 

—34,66—352t35 

4.817+jO.O 

-0.127 

-34 . 68  + j  52.29 

-34. 68-j52.29 

4, 7.51+j  0.0 

-0.162 

-34 . 69  + j  52  ,  27 

-34.69-j52,27 

4.733+j0.0 

-0.207 

-34 . 70  + j  52.25 

-34,7  0-j  5  2,25 

4 , 709  + j  0 . 0 
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previous  study,  some  values  of  g^3  caused  instability. 
Particularly,  for  gl3>~9,9,  pole  C  was  in  the  right  half 
plane,  The  greatest  sensitivity  was  experienced  for  values 
of  g13  near  g^3  =  -30.  Note  that  the  value  of  g  3  chosen 
via  Riccati  will  cause  stable  operation  and  is  also  removed 
from  the  value  of  greatest  sensitivity. 

3 .  Routh  Criterion 

The  Routh  criterion  was  applied  in  order  to  confirm 
the  range  of  feedback  gains  for  stable  operation  obtained 
from  the  root  locus  study.  The  Routh  array  is  set  up  below. 
The  C.E.,  given  in  Eq.  4,B,17,  is  used  with  the  understanding 
that  all  coefficients  are  absolute  values. 
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a21  +  a23  g2 


a33  +  g3 


('a23  gl  +  a21  S35 


l21  +  a23  S12  + 


a23gll"a21g13 


a3  3+S3 


-a 


2  3S1  +  a21S3 


The  conditions  for  stability  are: 
X)  §3  >  "a33 


( 4 , B . 14  ) 


since  a33  -  64.493,  then 

g3  >  -64,493  for  stable  operations. 
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2)  a21g3  >  a23Sl 


3) 
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If  g3  >_  -<64  the  g2  >  -201,72. 

The  boundary  limits  of  g2,  and  therefore  g  ,  are 
very  sensitive  to  g3< 

For  the  specific  control  chosen  in  Section  IVA,  the 
Routh  array  was  used  to  determine  the  stability  limit  for 
the  gains. 

From  Eq(4,B.14) 

s13  >  ~a33 

therefore  18.3  >  -'64,493 

From  Eq( 4 . B , 15 ) 

a21 

S11  >  S13 

therefore  1688  <  3689,581 

From  Eq(4.B.16) 

>  (a33Sll~f’a21a33) 

2  a23(a33+S13} 

therefore  119  >  -273.95 

It  was  seen  from  this  study  that  the  gains  were  well 
within  the  limits  of  stability.  Lower  gains  were  realitively 
more  stable. 
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V.  CONCLUSION 


It  has  been  shown  that  the  linear  regulator  design  can 
be  applied  to  the  problem  of  vertical  plane  motion  for  the 
captured  air  bubble  surface  effect  ship.  The  linear  pressure 
lift  only  model,  developed  from  the  simplified  equations  of 
motion,  was  shown  to  be  a  valid  model  which  could  be  used  for 
control  analysis  and  design.  An  optimal  control  was  found, 
based  upon  a  quadratic  cost  function,  using  empirically  de¬ 
rived  weighting  factors.  The  pressure  lift  only  linear 
system  showed  a  reduction  in  center  of  gravity  acceleration 

for  11  <ui  <100  radians/sec.  A  C.G.  acceleration  reduction 
e 

was  observed  in  both  the  linear  and  nonlinear  systems.  The 
shape  of  the  frequency  response  magnitude  curves  were  identi¬ 
cal.  A  13%  maximum  reduction  was  observed  at  to  =45  rad/sec. 

e 

where  peak  acceleration  occurs. 

The  control  gains  derived,  as  shown  in  Section  IV,  are 
within  the  airflow  rates  currently  available  for  the  XR-3. 

The  Riccati  equation  computer  program  on  Ref.  10  proved  to 
be  a  most  useful  tool,  which  facilitated  experimentation  of 
a  large  number  of  weighting  coefficients,  making  an  optimal 
control  design  more  feasible. 

Finally,  the  analysis  shows  that  the  system  is  stable 
over  the  range  for  the  control  gains  selected.  It  was  seen 
from  the  root  locus  study  that  the  limits  of  stability  narrowed, 
that  the  sensitivity  to  air  flow  rates  increased  and  C.G. 
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acceleration  decreased  with  higher  control  gains.  It  was 
seen  that  for  these  higher  gains  the  C,G,  acceleration  was 
not  reduced  and  the  control,  3  ,  increased  for  to  <17  rad/sec. 
This  was  due  to  the  fact  that  the  buoyancy  effects  are  pre¬ 
dominant  in  that  frequency  range  and  the  buoyancy  terms  are 
at  a  constant  steady  state  value  in  this  model. 

Though  an  optimal  control  was  found  based  on  a  particu- 

* 

lar  quadratic  cost  function,  many  other  possibilities  exist. 
Different  defined  weighting  factors  can  be  experimented  with, 
and  some  may  result  in  improved  C.G,  reduction.  The  limiting 
factor  in  such  a  design  is  equipment  limitation.  The  choice 
of  weighting  factors,  and  ultimately,  the  control,  are  even¬ 
tually  determined  by  that  fact. 
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VI.  RECOMMENDATIONS  FOR  FURTHER  STUDY 


The  following  list  of  recommendations  is  presented  for 
consideration  and  further  study. 

A.  A  more  realistic  system  presentation  could  be  accom¬ 
plished  if  pitch  coupling  was  combined  into  the 
linear  regulator  study.  This  would  facilitate  the 
speed  analysis  study  for  the  avoidance  of  wet  decks 
and  undesirable  pitch  accelerations. 

B.  Th'e  passive  control  technique  discussed  by  Boggio 

in  Ref.  6  should  be  examined  further.  The  introduc¬ 
tion  of  a  flexible  membrane  into  the  plenum  chamber 
may  reduce  vertical  accelerations.  This  technique 
has  been  examined  in  model  towing  tank  studies. 

C.  Since  the  mass  of  plenum  air  X3  cannot  be  directly 
measured,  an  observer  might  be  used  in  conjunction 
with  state  X3  or  perhaps  remodel  using  low  order 
model  analysis.  Perhaps  consideration  should  be 
given  to  a  design  using  only  XI  and  X2  as  states  to 
be  used  for  feedback, 

D.  The  linear  regulator  technique  design  approach  should 
be  compared  in  greater  detail  to  the  filter  design 
technique  proposed  by  Rohr  Marine  Inc,  in  Ref.  3.  A 
fixed  air  flow  rate,  based  upon  a  filtered  observed 
sea  state,  could  be  a  more  realizable  design  in  terms 
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r 

of  the  physical  equipment  limitations.  The  study 
of  fixed  air  flow,  of  Appendix  E,  shows  that  the 
filter  design  approach  may  offer  greater  control 
over  the  entire  encounter  frequency  range,  (see 
Fig,  E~l) 


i 

[ 
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APPENDIX  A 


The  simplified  non  linear  system  of  Ref.  1  is  shown  in 
Fig,  Al,  This  system,  with  the  buoyancy  term  removed,  served 
as  the  basis  for  study  in  this  thesis,  The  linearization 
and  resultant  linear  system  matrix  is  shown  for  the  simpli¬ 
fied  non  linear  system.  The  pressure  lift  only  model  lineari¬ 
zation  is  also  shown. 

LINEARIZATION  OF  THE  SIMPLIFIED  MODEL 

The  equations  of  motion  are  linearized  using  Taylor 
series  expansion  around  the  steady  state  operating  point. 
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Schematic  Diagram  of  Signal  Flow. 
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Writing  the  above  equations  in  state  equation  form; 
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The  characteristic  equation  was  found  from  the  system  matrix, 


DETCsI-A)  =  0 


to  be 


s  +  62s  +  2183s  +  789  =  0 


which  factors  into  $ 


(s+.37)(s+3Q.82  +  j34,8)(s+3Q,82-.j34,8) 


TAYLOR  SERIES  EXPANSION  OF  PRESSURE  LIFT  ONLY  VERTICAL 
EQUATIONS  OF  MOTION 
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then  incremental  acceleration  due  to  pressure  is 
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Then  the  incremental  volumetric  change  in  air  mass  is 
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A  series  of  system  coefficients,  k  ,  were  defined  for  ease 


of  writing  and  are  used  throughout  the  thesis.  They  are, 
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APPENDIX  B 


A  brief  step  weight  removal  study  was  performed  and 
compared  to  Ref.  1  as  an  additional  check  on  system  validity. 
The  IODE  remote  terminal  CP~CMS  program  on  the  IBM/ 360  was 
used  in  this  study.  A  sample  listing  is  shown  in  Appendix 
C. 
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APPENDIX  C 


The  following  are  the  listings  for  the  various  computer 
programs  used  in  the  preparation  of  this  thesis. 
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OPTIMAL  CGNTROL/ KALMAN  FILTER  PROGRAM 
PROBLEM  IOHNTI FICATI CM  -  HiAVS  CONTROL 


#***»«♦*»*******»»*»¥*« 4 * 4a 4*444  4************ 


THE  A  MATRIX 


0.0 

-2.267782000  C3 
-1.3CO2t>C0GC  54 

THE  e  MATRIX 


l.OOOJQOuJO  30 

8:8 


-i.i24309000  01 
-6.44S3Q0CQD  01 


0.0  C.Q 

THE  C  MATRIX 

-2.26778200C  03  0.0 


l.OOCQOOOOO  00 
-1.124800000  01 


**************«.**•**»  **********4  **«t«*****»*4 


***  CONTROL  OPTION  *** 


THE  R  MATRIX 
3. COO 00000 C— 02 
THE  C  MATRIX 


1.000C0003C 

05  0.0 

O.C 

O.C' 

1.3303 JCUOD 

03 

j  .0 

3.0 

C.O 

5.000030000-01 

STEAOY  STATE 

SOLUTION 

GAINS 

1 . 68839469  C 

03  -1. 191156  0  60 

02 

1.830578260  01 

**♦  CONTROL  OPTION  **• 


THE  R  MATRIX 

3.000030000-02 

- 

THE  0  MATRIX 

1 • 000033000  05 

3.0 

0.0 

O.C  , 

UCOCCC0000  03 

3.0 

0.0  1 

0.0 

S.CCCOOOOOD-Ol 

INITIAL  CCNOITI C  A S 

3.0 

0.0 

0.0 

0.0 

0.  ) 

1.9 

O.C 

0.0 

C.O 

****»• *4**4* 4**4 4 4*4^  4  44*4*4*444  4*44-»4  4**4444 


TIMe  2.CCC707C0C  Cl 

GA  INS 


TIME  *,  1.9CC00J32C  01 

GAINS 

1 .688395 Jo  0  33  -1. IS 1155710  12 

TIMc  *,  1.ECCOCOC40  01 

GAINS 

1. 688094690  03  -1. 191156360  32 

TIME  «,  1.7 CCCC0C6C  Cl 

GAINS 

1 • 638C9469C  03  -1.191156060  02 

TIME  «.  1.6CCCCCC60  Cl 

Ga  Ins 

1 .688094690  C3  -1.191156060  02 


0.0 

1.83057772D  01 
1.830578260  01 
1.830578260  01 


1. 833578260  01 


THc  INPOT  CATA  IS 
GRAPH  TITLC  AMO  SCALING  CATA 


ROOT  LCCLS  PLOT  FTP  "CCIFI'C 
PRESSUR '  ONLY  LINeAR  YCCEL 


XSCALP 

Y  SC  A  L  c 

IXUF 

I V  E IGF 

IWICF 

IFIGh 

IC-PTC 

20.0 

70.0 

3 

6 

7 

6 

I 

OROFF.  CF  TF?  Ch/ s  AC  TCR  I  ST  IC  CCLATICN 
3 

CONSTANT  PART  OF  T I- 2  COEFFICIENTS  IN  CSSC^MOtNG  CPOrR 
O.IOCCCO  01  C.O  0.2267EC  CA  C.5T010C  CO 

VARIASLF  FART  OF  TF'  COEFFICIENTS  IN  CESC^NC  INC  CPDCR 

0.0  0. 10 AS 70  02  c.a  c.c 

INITIAL  VALLc  Cl-  71-2  WFIAeiS 
G.500CGE  01 

NUMBER  OF  CFCA'JFS  TO  3C  SPANflSC 

l 

THE  SYSTEM  POL'S  A  Sit 

»**  I  F  S  UU'-^TSTI  r*»  WARNING U ITF  FIX’)  2QA0R  2  (ICR  *  66) 

REAL  PART 

-0.25 

C.C 

REAL  PAST 

O.;)  0.0 

IMAGINARY  PART 


3.1 2  570F-C3 
imaginary  part 


J.12570S-C3 


0. A762 lc  C 2  -*.A7<21e  C2 

THE  SYSTEM  2 FRO S  ARC 


0.0  C.C 

ALL  POINTS  HAW"  Tr  -  SAME  CCCFC IN  AT5S . 
THIS  PLOT  WILL  NOT  OF  CLTFLE. 
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- SPECIFICATIONS - 


VARIABLES  £  INITIAL  CONEITIONS: 

XI  *  0.0 
X2  *  C.C 
X3  »  0.0 
T  *  C.C 

CONSTANTS: 

A21  -  2267.782CC0 
A 23  *  U.248C0CQ0 
A 2  1  »  2116.265CCC 
A33  *  1O.497C00CC 
621  >  2267.782000 
621  *  21 1 6  .2  65  COO 
832  a  l.CCCOCOCOO 
G  »  32.174C0C00 
A  »  C.1CCCCCCCCGD  00 
HE  »  45 .5CCCGCCC 
K6  a  6.14422C0C0 
K1 1  »  -1666.C94700 
K12  >  119. 1160C00 
K13  «  -16. 2C6CC000 

SPECIAL  FUACTICNS: 

LM  «  (2.*3.141S9265*Gl /WE 
r  *  CSIM62.63I9/L("J)/(6Z.8319/LWJ 
SEA  a  A*SINIWE*T> 

U2  »  K11*X1+X12*X2+K13*X3 
V  a  -A21*Xl-A23*X3-S2l*F*SEA 

DERIVATIVES: 

C*X1  /C(T  J  »  « 

X2 

01X2  /OIT  )  »  » 

-A2i*Xl— A23*X 3—821* F* SEA 
0CX3  AC(T  1  a  * 

-A3l*K6*X l-A33*K6*  X3-B31 *K6*F*SE A+U2 


OUTPUTS: 

TITLE:  CAB  SES 
TABULATE:  T 
AT  INTERVAL 
NC  PL CT TING 


INCREASEC  AIR  FLOW  RAT 
Y 

0.5CQ0CO00 000-02 


ENC  CALCULATION  WHEN  T  „GE.  1.50000 
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SPECIFICATION 


VARIABLES  £  INITIAL  CONDITIONS; 

XI  =  0.500000 
X2  =  0.0 
X3  =  "44.5300 
T  =  0.0 

SPECIAL  FUNCTIONS: 

Y  =  2030.26*X1J.0.09*X3 

DERIVATIVES: 

D(X1  /D(T  )  =  = 

X2 

D(X2  /D(T  )  =  = 

-2182 . 68*X1-10 . 76*X3 
D(X3  /D(T  0  =  = 

-12504. 36*X1-62.0*X3 

OUTPUTS : 

TITLE:  SES 

TABULATE:  T  XI  X2  X3  Y 

AT  INTERVAL  0.300000D-01 
PLOT:  X2 

AGAINST:  T  AT  INTERVAL  0. 10000 0D- 01 

END  CALCULATION  WHEN  T  , GE .  0.600000 
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APPENDIX  D 

The  characteristic  equation  for  the  pressure  lift  only 
model  under  control  conditions  was  derived  from  the  system 
matrix  with  the  feedback  control,  gains  added  to  the  state 
.  The  result  was  compared  with  the  characteristic  equation 


The  system  matrix  is 


p 

1 

0 

> 

11 

pi 

N> 

0 

a23 

_a 31 

0 

^33 

and  with  feedback  gains ; 


0  10 

A  =  a21  0  d23 

a3l+gll  +S12  a33+S13 

The  characteristic  equation  with  gains  becomes 

S  "(a33+g13)s  "(a23S12+a21)s 
"a23a3l"a23Sll 
+a2la33+a21S13 


It  should  be  noted  that  the  above  equation  is  identical  to 
that  derived  in  Ref,  1  except  that  the  gains  are  shown  in 
(4.B.4) 

a21  =  -2267.782 
2  =  ,248 

a31  =  -13002.8 
a33  =  -64,49? 
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From  the  Riccati  equation,  it  was  seen  that  the  gains  had  the 
following  sign  values ; 


"gll  *  +S12  1  "s13 

Substitution  of  the  appropriate  signs  yields 

s3+(a33+Sl3)s2+(a23§l2+a21)s 

+a21a33+a21S13 

"'a23a3l''a23Sll 

Substituting  in  the  system  value  parameters,  the  character¬ 
istic  equation  becomes: 

s3  +  (64,49  3+gn  3)s2+(  2  267.7  8  2+11.24  8g12)s-11.24  8g11 


APPENDIX  E 


FIXED  AIR  FLOW  RATE  STUDY 

A  brief  study  was  conducted  to  examine  the  effects  of  kc 

D 

on  the  overall  frequency  response  and  compared  with  the  op¬ 
timal  control  design  results.  Earlier  work  indicated  that 
the  system  response  was  sensitive  air  flow  rate.  Therefore 
three  values  of  kg  were  tested  on  the  pressure  lift  only 
non  linear  system  with  the  resultant  response  curves  plotted. 
Values  for  k~  are : 

b 

k-'  =6.144  ,  Fig.  El 

kg  =  12.0  ,  Fig.  E2 

kg  =  19 , 0  m  Fig,  E3 

Note  that  the  high  acceleration  midband  frequencies  are 

attenuated  as  k_  increases.  The  response  curve  flattens 
o 

out  and  midband  frequencies  are  reduced  by  20%.  lis  result 
would  tend  to  confirm  the  contention,  of  Ref.  3,  that  fil¬ 
tered  fixed  air  flow  rates  method  of  control  may  have 
advantages.  The  time  response  for  sinusoidal  input  remains 
linear  throughout  the  range  of  interest.  The  time  versus 
acceleration  plot  shown  is  a  time  history  for  an  air  flow 
rate  of  Qjq  =  60  at  an  encounter  frequency  of  wg=45  rad/sec. 
Appendix  C  presents  the  DSL  Bode  program  listing.  Finally, 
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NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA  F/G  13/10 

SURFACE  EFFECT  SHIP  HEAVE  CONTROL  USIN6  A  LINEAR  REGULATOR  DESl— ETCCU) 
MAR  80  D  H  EVERETT 


the  system  peak  frequency  magnitude  plot  for  the  conditions 
of  no  control,  state  variable  feedback  control  and  fixed  air 
flow  rate  control  is  shown  in  Fig,  E8  for  QjQ=60. 
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.00  -16.00  -8.00  0.00  8. 00  16.00  24.00  32.00  40.00  48.00 


CAB  SES  K S  =  S . 14422 

BODE  PLOT,  MAG  (DB)  VS.  LOG  W 


XSCflLE=  0.40  UNITS/INCH 

YSCfllE=  8.00  UNITS/INCH 

Figure  El 


96 


CRB  SES  K6= 1 2 . 0 
gBODE  PLOT ,  MRG(DB)  VS.  LOG  N 


XSCflLE=  0.40  UNITS/INCH 

TSCflLE=  8.00  UNITS/INCH 

Figure  E2 


,-32.00  -24.00  -16.00  -8.00  0,00  8,00  16.00  24.00  32.00  40.00 


CRB  SES  K 6  = 1 9 . 0 

BODE  PLOT,  MRG  (DB)  VS.  LOG  W 


XSCRI_E  =  0.40  UNITS/INCH 

YSCRLE  =  8.00  UNITS/INCH 

Figure  E3 
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XSCRLE  = 
TSCRLE  = 


0.40  UNITS/INCH 

8.00  UNITS/INCH 


Fiugre  E6 
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